Asymptotically exact dispersion relations for collective modes in a confined charged 

Fermi liquid 
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Using general local conservations laws we derive dispersion relations for edge modes in a slab of 
electron liquid confined by a symmetric potential. The dispersion relations are exact up to X 2 q 2 , 
where q is a wave vector and A is an effective screening length. For a harmonic external potential the 
dispersion relations are expressed in terms of the exact static pressure and dynamic shear modulus 
of a homogeneous liquid with the density taken at the slab core. We also derive a simple expression 
for the frequency shift of the dipole (Kohn) modes in nearly parabolic quantum dots in a magnetic 
field. 
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The examples of importance of a collective motion of 
particles range from such old problems as energy losses by 
fast particles in solids to modern studies of a collective 
response in semiconductor nanostructures and interac- 
tion of clusters and molecules with intense laser radiation. 
However, only a few exact features of collective dynam- 
ics in spatially inhomogeneous many-body systems have 
been found up to now. The most known result is the Gen- 
eralized Kohn Theorem (GKT) [1] or the Harmonic Po- 
tential Theorem in a more general formulation [2]. This 
theorem has been extensively used for construction of ap- 
proximations for exchange-correlation (xc) potential in 
Time Dependent Density-Functional Theory (TDDFT) 
(see e.g. Ref. 3 and references therein). 

In this paper we present another asymptotically exact 
result which, in some sense, can be viewed as an exten- 
sion of GKT. We derive the exact dispersion relations 
for the dipole and monopole plasma modes in a slab of 
a charged Fermi liquid confined by a symmetric external 
potential. This result opens a possibility for a direct ex- 
perimental determination of the dynamic shear modulus 
of an electron liquid. It can also help to control the ac- 
curacy of approximations used to describe dynamics of 
inhomogeneous systems. It has also direct implications 
for the theory of surface plasmons in simple metals. In 
addition we derive a few asymptotically exact relations 
for dipole modes in quantum dots (QD) and wires in a 
homogeneous magnetic field. 

The collective motion of an arbitrary interacting sys- 
tem obeys the exact local conservation laws: 



d t n + d^riVf! = 0, 
mnDtVft + en(v x H) M 



(1) 
0, (2) 



where D t = dt + v v d v , n is the density of particles, 
v = j/n is the velocity, U = J7 cx t + Un is a sum of the 
external and Hartree potentials, H is the external static 
magnetic field and P^ v is the exact (generally unknown) 
stress tensor. 

Introducing the displacement vector u as = dtu^, 



we derive from Eq. (2) the linearized equation of motion 

mn <9 2 u M + en (d t u x H) M 

+ d v SP^ + n d IM SU + Snd ll Uo = Q, (3) 

where Sn = —d^n^u^, 8P^ V and SU are deviations of the 
density, the stress tensor and the potential from their 
equilibrium values no, P°2 an d Uq respectively. Eigen- 
frequencies of Eq. (3) define the exact excitation energies 
of a system. 

Let us first consider the dipole modes of a finite system 
confined in D dimensions by the potential t/ cx t m the 
presence of a homogeneous magnetic field H. Integration 
of Eq. (3) leeds to the following exact equation for the 
dipole moment per particle d — J un dx 



muj d,. 



iwe(d x H) M - i J noUud^duU^dr. (4) 



This equation, commonly called the dynamical force sum 
rule [4] , shows the exact decoupling of the center-of-mass 
and relative motions for an external potential of the form 
Ucxt = \ m ^u,v x i^ x v [1; 2]. In this case the dipole (Kohn) 
modes correspond to the rigid motion with the displace- 
ment u = d = const being the solution to the Newton 
equation [2] 



iuje(d x H) A 



(5) 



In practice one often deals with the quasiparabolic QD, 
where the deviation AU cx t of the confining potential from 
the parabolic one is small (see Ref. 5 and references 
therein). The perturbative treatment of Eq. (4) shows 
that to the first order in At/ cx t the dipole moment still 
satisfies Eq. (5) but with the renormalized dynamical ma- 



trix nl„ 



A17 2 Ily , where 



= N 



no 



dx^dx^ 



AUcvtdr 



(6) 



The corrections to the Kohn-mode frequency has been ex- 
tensively studied experimentally and theoretically (using 
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approximate many-body methods) [5]. The importance 
of Eq. (6) in this context is that it exactly accounts for the 
many body effects at the linear in A[/ ex t level. Therefore 
it can be helpful for studying correlation phenomena such 
as Wigner crystallization in a few-electron QD (see e.g. 
Ref. 6) by means of far-infrared photoabsorption. The 
transition to the correlated state, which can be tuned 
by a magnetic field, changes the density distribution no 
and consequently affects the frequency shift according to 
Eq. (6). 

Another useful relation for a dipole mode can be ob- 
tained for a system of N charged particles confined by 
a D-dimensional potential (D= 3, 2, 1) with a spherical 
(dot), cylindrical (wire) or a slab geometry (we consider 
the case H = 0). We assume that n (r) acquires a con- 
stant value n inside a system and vanishes at infinity. 
The length Ao of the density variation at the boundary 
is assumed to be smaller than the system size o, which 
is defined in such a way that the rectangular density dis- 
tribution n0(a — r) provides the correct number of parti- 
cles N. Noting that the density fluctuation Sn(r) for the 
dipole mode is localized at the length scale As <C a near 
the system edges, we obtain from Eq. (4) the following 
result 

where Sn(r) = Sn(r)/ 5n(r)dr is the normalized ra- 
dial part of the density fluctuation and AJ7 ext = U cxt — 
^har> where the harmonic potential C/har = mui 2 r 2 /2D 
corresponds to the same equilibrium density in the inter- 
nal region (ui 2 = Aire 2 n/m). It is noteworthy that A[/ ox t 
can be arbitrary lage in the edge region as the smallncss 
of the right hand side in Eq. (7) is controlled by As /a. 

Equation (7) describes the the frequency shift relative 
to the classical value uj p /VD. The shift is either positive 
or negative depending on whether [/ ox t is "harder" or 
"softer" than Uhm- It obviously vanishes for A£/ ext = 
in accordance with GKT. Since 8n{z) has an imaginary 
part, the collective mode decays due to Landau damping. 

For a wire or a slab geometry Eq. (7) gives the fre- 
quency at zero wave vector q. Now we consider a general 
case q^O, which requires a local treatment of the basic 
Eq. (3). 

Consider a slab with a given sheet density N s con- 
fined in z-dircction by a symmetric potential U ex t{z). 
The equilibrium stress tensor can be written as P^ = 
Po(z)S fJ , u + TT < ^,(z), where P = |TrP°3 is the exact pres- 
sure and 7r^ is a traceless tensor which describes quan- 
tum effects. The symmetry requires that P^J — for 
j = x, y and Pg = P ± = P + tt , P°% = Pfl = P\\ = 
Po - ttq/2. 

The equilibrium density no{z) satisfies the balance 
equation (the static version of Eq. (2)) 

d z P c J + n d z U = 0, (8) 



where Uq = J7 cxt + Uu and Uu is related to n by the 
Poisson equation 

8%U U = -47re 2 n . (9) 

Combining the relation Sn = —d^n u^ and Eq. (3) and 
considering the plane wave solution e ^ J ("*~ t i r ii) we arrive 
at equation for the density fluctuation <5n q ^(z) (indexes 
q, uj are suppressed below) 

muj 2 Sn + d 2 z 5P zz - qkqjSPkj + 2qjd z 5P zj 

+d z 5nd z U + d z n d z 5U - q 2 n 5U = 0, (10) 

where SU satisfies the equation 

(d 2 - q 2 ) SU(z) = -4ire 2 6n{z). (11) 

Up to this point all transformations are exact. From now 
on we consider the equilibrium density no(z) that satisfies 
the assumption used in the derivation of Eq. (7) [7]. 
For a slab geometry the width 2a is defined as follows 

/oo 
n (z)dz = 2an = N s . (12) 
-oo 

Let us choose the origin at the slab "edge" so that the 
center of the slab is at z = —a. Equation (12) im- 
plies that the density n (z) at z > —a can be repre- 
sented in the form no(z) = h9(—z) + Ano(z), where 
J^° a Ano(z)dz = 0. The function Ano(z) is localized at 
the length scale A . 

Our assumption allows to determine the asymptotic 
form of the stress tensor in the slab interior where 
n (z) rts n and the gradients of the density vanish. Here 
7r^ w and the equilibrium stress tensor reduces to the 
pressure P^ w PqS^ (all quantities marked with a bar 
refer to the slab interior). In the same region the lo- 
cal noncquilibrium correction to the stress tensor SP^ is 
defined by the asymptotic expression 

5Pp V = -KuS^du - [i u {dy,u v + <9„?v - ^S^duj 

(13) 

Within the limits of Landau Fermi liquid theory the bulk 
(K u ) and the shear (p,^) moduli are real and can be ex- 
pressed in terms of Landau parameters [9, 10]. In general 
K u and p, u have imaginary parts which are related to in- 
ternal viscosities [9]. 

Now we focus on two collective modes at the edges 
of the slab. These are the lowest antisymmetric (dipole) 
and symmetric (monopole) plasma oscillations. At q — > 
the dipole mode corresponds to the oscillations across the 
slab (which is the Kohn mode in the case of a parabolic 
confinement), whereas the monopole mode is a periodic 
in g-direction compression and expansion of the slab (2D 
plasmon). At 2qa > 1 the both modes merge into two 
surface plasmons. 
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We derive the dispersion relations for these modes with 
an accuracy A 2 g 2 for arbitrary values of the parameter 
aq, where A = max{A, Xg}. Let us define an intermediate 
length scale a > I > A and integrate Eq. (10) over the 
region (— l,oo). The result of the integration is 



_p_ 

2 

2 r 



(1±. 



-2qa\ 



1 — q(z8ri) — 7r(z 2 8n) 



Q_ 

m 



2jlu 
n 



+ 



(Sn) 



D 



+ O(AV) 



(14) 



where (...) = ...dz and 8n(z) = Sn(z)/(Sn). The 
upper(lower) sign corresponds to the dipolc(monopolc) 
mode and the function D is defined as follows 

D = 2ne 2 {n(0{-z)z 2 5n) 

+ ((\z-z'\An (z)SnJz')))Te~ 2qa (zAn )} (15) 

Since both Auq and Sn vanish at — z ~ A < I, all integrals 
in Eqs. (14) and (15) are independent of the lower limit. 
Equation (14) is, in fact, an implicit dispersion relation. 

To simplify integrals (zSn) and (z 2 Sn) in Eqs. (14), 
(15), we multiply Eq. (10) by z and z 2 respectively and 
integrate over z from — I to oo. After calculations we 
arrive at the following equations 



m(uj 2 
m(uj 2 /2 



uj 2 )(zSn) = (Snd z AU cxt ) 

4ne 2 q (1 T e- 2qa ) (zAn ) + 0(A 2 a 2 ),(16) 



- u 2 p )(z 2 Sn) = -(SP ± )/(Sn) 
+ (zfrnd z AU cxt )-D + 0(\q), 



(17) 



Using Eqs. (14)-(17) and introducing the notation 

f± = l± e~ 2<l \ 
we obtain the dispersion relations in the final form 

4 = y/± + ^(Snd z AU cxt ) + ^(A T +B T ) (18) 
with the following coefficients in the second-order term 
"\/2/± 



(Snd z AU cxt ) - (zSnd z AU cxt ), (19) 



= 2p u /n+((5P [ \+6P ± ))/(5n) 
+ 2ire 2 (2(zAn )f ± - (z 2 5n)fe) 



(20) 



Let us analyze Eqs. (18)-(20) in more detail. 

If U cx t differs from the harmonic potential the main 
contribution to the dispersion comes from the first two 
terms in the right hand side of Eq. (18). In the limit 
2qa <C 1 the lower mode becomes a 2D plasmon with the 
universal dispersion io\ = 2ire 2 N s q/m. On the contrary, 
the frequency of the higher (dipole) mode 



contains a non-universal shift which is proportional to 
the inverse width of the slab. At q — Eq. (21) recovers 
the result of Eq. (7) for D=l. However a close inspection 
of Eqs. (14)-(17) shows that, for a dipole mode, they are 
consistent only for q > A/a 2 , which seemingly makes the 
limit q — > inaccessible. Yet in the region 1 > qa > 
A/a, where Eqs. (14)-(17) are applicable, the dispersion 
Eq. (21) is a linear function of q which is shifted by a 
constant as given in Eq. (7). Thus, despite Eqs. (14)- 
(17) are ill defined at q — > 0, the resulting dispersion 
relations, Eq. (18), are correct up to q = 0. 

In the opposite limit 2ga 3> 1 (/± — ► 1) both modes 
merge into antisymmetric and symmetric combinations 
of surface plasmons at two slab surfaces with the same 
frequency w T = u) S 



= w p/ 2 + q{5nd z AU cxt )/m. 



(22) 



In a particular case of the jellium edge d z AU cxt = 
—muo 2 z6{z) and we recover the result 



= w p/ 2 - qu p 



f 

Jo 



zSndz, 



(23) 



obtained in Ref. 11 using the dynamical force sum rule. 

The intriguing feature of Eqs. (18)-(20) is a very special 
role of the harmonic potential. Both the linear term in 
Eq. (18) and one of the second order coefficients A T , 
Eq. (19), vanish if AU cxt (z) = [12]. This is obviously 
related to GKT [1]. The existence of the exact GKT 
mode with oj = uj p at q = can be easily observed via 
direct substitution of the rigid-shift solution: u^(r,t) — 
S^ z u z (t), and 

5n{r,t) = -u z 8 z n , <5P^(r,t) = -u z d z P$ (24) 

into Eqs. (3) and (11) at q = with i7 ex t = Uhm- In- 
terestingly, the solution Eq. (24) remains asymptotically 
valid at finite q for both modes under consideration. Sub- 
stituting Eq. (24) with u z = u z (t)e iqr ii into Eqs. (3) and 
taking into account the Poisson equation for SU we ob- 
tain the equation 

d^u z {t) + 27re 2 n (l ± e" 2<?a ) u z {t) + O(qz) = 0, (25) 

which shows the validity of the rigid-motion solution in 
the region z <C l/q near the edges of the slab. This 
observation allows to simplify the remaining coefficient 
B T [Eq. (20)] in Eq. (18). Using Eq. (24) we find that 
the integrals of the stress tensor fluctuations (the second 
term in Eq. (20)) are reduced to the equilibrium pressure 
per particle in the internal region 



(SP±)/(6n) = (SPu)/(Sn) = P /n, 



(26) 



whereas the integral (z 2 Sn), which enters the last term 
in Eq. (20), is transformed to the form 



-qaul + (Snd z AU cxt )/am (21) 



2ire 2 n(z 2 Sn) = Aire 2 (zAn ) . 



(27) 



4 



The right hand side of Eq. (27) is, in fact, the potential 
of a double layer which is formed by the density distri- 
bution An a (z). On the other hand, the integration of 
the balance equation Eq. (8) leads to the following exact 
relation of the double layer potential to the pressure per 
particle inside the slab 

4ire 2 (zAno) = P /n. (28) 

Substituting Eqs. (26)-(28) into Eqs. (20) and (18) we 
obtain explicit dispersion relations for collective modes in 
a parabolically confined system 

u T = -*-f± + —f± + — — q • 29) 

T 2 \mn ran ) 

In the surface plasmon regime (2qa > 1, oj t = ojs) this 
equation simplifies further as 

lo 2 = uj 2 p /2 + 2q 2 (P + ft,) /mn = uj 2 p /2 + v 2 s q 2 . (30) 

This equation relates the dispersion of plasma modes to 
the exact pressure Pq and shear modulus p, u of a ho- 
mogeneous electron liquid at a given temperature T and 
density n. Equation (29) shows the absence of Landau 
damping up to 0(X 2 q 2 ). The attenuation is due to the 
internal viscosity (the imaginary part of p, u ) which is re- 
lated to the multi-pair excitations. From exact results 
Eqs. (29) or (30) we can easily obtain the dispersion in 
any particular approximation. For example, in Random 
Phase Approximation (RPA), which assumes Hartree ap- 
proximation at the static level, we have p, u = Pq [10] 
and, therefore, v\ = iP /mn (i.e. |uf. at T = 0). This 
is times 4/3 of the corresponding coefficient in the RPA 
bulk plasmon dispersion. 

As the exact pressure Pq (n) is presently available from 
Monte-Carlo calculations, Eq. (29) opens a direct access 
to the complex dynamic shear modulus p, w in the range 
< uj < ujp/\/2. To determine p,^ one has to mea- 
sure the dispersion of plasma modes in a wide parabolic 
quantum wells using e.g. the grating technique [13]. To 
access the quadratic part of the dispersion, which is gov- 
erned by the shear stress, one needs thinner grating and, 
possibly, wider wells than those used in Rcf. 13. Such 
measurements would allow a direct comparison to recent 
many-body calculations [14]. 

As any exact result, Eq. (29) should help to control 
the consistency of approximate methods such as differ- 
ent hydrodynamical approaches [10, 15, 16] or TDDFT 
schemes with approximate xc kernels / xc [3] . 

Let us consider hydrodynamics first. Since the dis- 
persion Eq. (29) is controlled by the shear stress, any 
hydrodynamics with a diagonal stress tensor SP^ ~ <5^„ 
[15, 16] does not reproduce the correct dispersion co- 
efficient. In the adiabatic hydrodynamics [15] also the 
dispersion of bulk plasmon is wrong. This problem is re- 
moved in the theory of Rcf. 16, which is however valid 



only for ID motion, where the tensor structure of SP^ 
is irrelevant. The generalized hydrodynamics of Ref. 10 
gives the exact result since it reproduces the correct 
structure of 5P^ U . 

Similar arguments apply to different approximations 
for f xc in TDDFT. The kinetic part of /x, which is equal 
to the pressure of Kohn-Sham particles, is reproduced 
correctly. The xc part is, however, / xc -dependent. For 
example, the adiabatic approximation gives no additional 
xc contribution to the shear stress which is completely a 
non-adiabatic effect. The consistent result can be ob- 
tained using the approximate xc kernel of Ref. 8. The 
correct asymptotic form of the stress tensor Eq. (13) was, 
in fact, one of the requirements in the derivation of this 
approximation. The dispersion of surface plasmons is an 
example of physical situation where this asymptotic re- 
quirement leads to an observable effect. 
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